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I. INTRODUCTION 


There are quite a lot of empirically formulated equations to 
calculate the water wave characteristics at the point 
breaking wave. These equations are called the breaker 
index. There arebreaker height index to calculate breaker 
height,breaker depth index to calculate breaker depth with 
an input of breaker height, andbreaker length index to 
calculate breaker length with the inputs of breaker height 


and breaker depth. 


Some researchers who put forward equation for the 
breaker height index include Le Mehaut and Koh (1967), 
Komar and Gaughan (1972), Sunamura and Horikawa 
(1974), Singamsetti and Wind (1980), Ogawa and Sutto 
(1984), Larson and Kraus (1989a), Smith and Kraus 
(1990), Gourlay (1992), and Pittikon and Shibayama 


(2000). 


Meanwhile, studies dealing with breaker depth index are 
conducted by Mc Cowan (1984), Galvin (1969), Collins 
and Weir (1969), Goda (1970), Weggel (1972), Sunamura 
(1980), Singamsetti and Wind (1980), Seyama and Kimura 
(1988), Larson and Kraus (1989b), Smith and Kraus 


(1990). 
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Abstract— In this study, equations to calculate breaker height, breaker 
length, and breaker depth were formulated. The analytical formulations 
were based on the velocity potential of the solution of the Laplace equation. 
The substitution of velocity potential to weighted kinematic free surface 
boundary condition yielded equation for wave amplitude function. The 
dispersion equation was obtained by substituting the velocity potential and 
wave amplitude function into the weighted momentum Euler equation. The 
relation between breaker height and breaker wavenumber is extracted from 
the wave amplitude function. The breaker wavenumber equation was 
obtained by substituting this relation into the dispersion equation. 
Furthermore, the breaker depth equation obtained was put to the 
wavenumber conservation equation with the variables of breaker height 
and breaker wavenumber. The breaker wave amplitude equation, the 
breaker wavenumber equation, and the breaker depth equation are simple 
equations that are very easy to use. 


In terms of breaker length index, researchers proposing the 
studies include Miche (1944), Battjes and Jansen (1978), 
Ostendorf and Madsen (1979), Battjes and Stive (1985), 
Kamphuis (1991), Rattanapitikon and Shibayama (2000) 
and Rattanapitikon et al. (2003). 


This study formulated equations for 
characteristics at the point breaking wave based on 
hydrodynamic equations including the equations for 
breaker height, breaker length, and breaker depth. Using 
this method, a relationship among the three is elaborated. 


water wave 


In this research, the wave amplitude function was 
formulated using the weighted kinematic free surface 
boundary condition and the weighted Euler’s momentum 
equation. Hence,there is a weighting coefficient in the 
equations of breaker wavenumber, breaker depth, and 
shoaling in the dispersion equation. The weighting 
coefficient is obtained by calibrating the critical wave 
steepness with the critical wave steepness criteria from 
Michell (1893), where the critical wave height used is the 
maximum wave height from Wiegell (1949; 1964). In 
addition to calibrating the critical wave steepness, the 
weighting coefficient is also obtained by calibrating the 
breaker height with the breaker height from the breaker 
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height index equation. The breaker depth obtained is also 
calibrated to the breaker depth from the breaker depth 
index. 


The breaker length generated from the breaker length 
index equations has a very wide spread and is much 
different from the critical wave steepness either from 
Michell (1893) or Toffoli et al (2010). Therefore, breaker 
length calibration was not carried out in this study. Only 
comparisons between the breaker length model and the 
breaker length from the Miche’s (1944) breaker formula 
were made. 


Il. VELOCITY POTENTIAL 


The complete velocity potential obtained from the solution 
of the Laplace equation with the separation variable 
method (Dean, 1991) is: 


o(x,z,t) = G(coskx + sin kx)coshk(h + z) sino t 


k is the wavenumber, in general, k and G are wave 
constants,o is the angular frequency, xis the horizontal 
axis at the still water levels, while z is the vertical axis, 
where z = Oat the still water levels. 


In this equation, there are two components of velocity 
potential. They are thecos kx and sin kxcomponents. The 
constants G and kmust have the same value on the two 
components of velocity potential. Therefore, the two 
constants are determined at a value of kxwhere the two 
sinusoidal functions are the same. The point is called the 
characteristic point. The characteristic point on the 
horizontal xaxis is the point where the value of cos kx is 
equal to the value of sin kx. Likewise, the wave amplitude 
which is also a wave constant is determined at the 
characteristic point. 


At this characteristic point, analysis can be done using 
only one component of the velocity potential, which in 
this study is the cos kx . 


@(x,z,t) = G cos kx coshk(h + z) sin o t....(2) 


The use of one component implies that the value of G in 
(2) is the sum of the components cos kx and sin kx, or that 
G in (2) has double values. 


The water particle velocity in the horizontal x direction is, 


d 
u(x,z,t) = — 7 = Gk sin kx coshk(h + z) sino t 


(3) 
The water particle velocity in the vertical z direction is, 
do 
w(x,z,t) = — = —Gk cos kx sinhk(h+ z) sino t 
Z 
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(4) 

In the velocity potential equation, there are two 
conservation equations (Hutahaean, 2021), which are the 
energy-conservation equation and the wavenumber 


conservation equation. 


a. Energy Conservation Equation 





GÆ +2k=0 un (5) 

dx dx 
b. Wavenumber Conservation Equation 
ak(h+5) 

= 0.....(6) 


In whichAis the wave amplitude. 


MI. WAVE AMPLITUDE FUNCTION 


The wave amplitude function was formulated using the 
Weighted Kinematic Free Surface Boundary Condition 
(KFSBC) from Hutahaean (2021): 


dn = _ dn 
Van = Wn Un a) 


y is the weighted coefficient that is positive and greater 
than 1, 7 = n(x,t) is the water wave surface elevation 
equation, w, is the particle velocity in the vertical 
direction on the water surface, while Un is the water 
particle velocity in the horizontal-x direction. 
Substitutions (3) and (4) on (7), 

dn 


Yar = —Gk cos kx sinhk(h + n) sino t 


d 
—Gk sin kx coshk(h +n) sino t 7 


Forcoskx # 0, 





dn sinkx dn 
Yar = —Gkcoshk(h +n) (tanhk(h +n)+ ) 


cos kx sino t 


cos kx ax. 


For a periodic function, 





Gkcoshk(h + (t hk(h + a 
coshk(h + n) |tanhk(h +n) zos kx Ax 


constant 
This equation is a wave constant, whose value is calculated 


at the characteristic point, the point where coskx = 
sin kx, 


a 
r= = —Gkcoshk(h + n) (tanhk(h +n) 


dn , 
+ A cos kx sino t 
dx 


This equation is integrated to time-t, 
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Gk 
n(x,t) = mao +n) (tanhk(h +n) 
dn 
+ A cos kx coso t 
dx 
The wave amplitude function defined is, 
a cote: (t hk(h + +20) 
= yos ( n) | tanhk ( n) ae 


Given that G is doubled, the right-hand side must be 
divided by 2, 


eer eee Y(t hk(h + +5) 
aap n) (tan DF 


The water level equation becomes, 
n(x,t)=Acoskxcosot .... (8) 


The wave amplitude, which is also a wave constant, must 
be constant and have the same value as the cos øt and 
sin o tfunctions. Thus, the value of the wave amplitude is 
determined at the characteristic points x and t where, 
cos kx = sin kx and cos ot = sinot. 


At this characteristic point, 


A an kA 
n = -and— = —— 
2 dx 2 


The equation for the wave amplitude function becomes, 


a= casa n+5)(¢ hk n+5) S 
Caa (eee 3) 


Given that the wave amplitude A is constant, thus: 


A 
tanhk (a + 5) = constant 
and 

A 
coshk (a + 5) = constant 
Thus, 

A 

k (n + 5) = constant 


Defined as, 
A 
k (h + $) = 6r...(9) 


Ois a coefficient whose value will be determined by 
calibrating the breaker depth to the breaker depth from the 
breaker depth index. It is defined as: 


Cy = tanh @7....(10) 
B=cosh@x dan £, = sinh @z....(11) 


Wave amplitude function becomes, 


= BOE (cn z =) ve (12) 


2yo 2 
This equation can be written as an equation for G: 
2yoA 


G= Bic(cn 4) ea cute (13) 
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IV. DISPERSION EQUATION 


The dispersion equation is formulated using the weighted 
Euler surface momentum equation (Hutahaean, 2021): 
1d 


du, = an 
poao ti (yuyu, + Wwa) = —9 ae 4) 


Substituting u from (3) and w from (4), 7 of (8) with A 
from (12), and the obtained equation is worked out at the 
characteristic point. 





> 1 og . gk kA 
ya + (B-B) 5 Gk? => (cn - >) 


Then, G is substituted by (13) and the equation is 
sic kA 
multiplied by(cn = ~), 


y7o (cn — 5 + (y —cf)yokA = = (cn — a5 


This equation is a dispersion equation, with the wave 
amplitude A as the variable whose value is known. 


To get the appropriate y weighted coefficient, Michell’ s 
(1893) critical wave steepness criteria are used: 


H 0142 ...... (16) 


-= 
The wave amplitude at (15), the maximum wave height at 
a wave period of Wiegel (1949; 1964) is used: 


za ous (17) 


~ 15.62 





Assuming a sinusoidal wave, then A = Z, In (17), g is the 


gravitation (m/sec?), T wave period (sec). The results of 
calculations with the value of the weighting coefficient 
y = 1.211 and 0 = 0.625 are presented in Table (1), with 
a critical wave steepness referring to Michell’s (1893) 
criteria. 


Table.1: Wavelength and Critical Wave Steepness in Deep 

















Water. 
T Ho Lo Ho 
(sec) (m) (m) Lo 
6 1.451 10.222 0.142 
7 1.975 13.913 0.142 
8 2.58 18.172 0.142 
9 3.265 22.999 0.142 





10 4.031 28.393 0.142 



































11 4.878 34.356 0.142 
12 5.805 40.887 0.142 
13 6.812 47.985 0.142 
14 7.901 55.651 0.142 
15 9.07 63.885 0.142 
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The wavelength in (15) is the wavelength at deep water ho 
which can be calculated by (9) as follows, 
On 


A 
ho = = Aol) 


Table.2: Deep water depth hy and “2 
0 



































T Ho ho ho 
(sec) | m) (m) Lo 
6 1.451 2.831 0.277 
7 1.975 3.854 0.277 
8 2.58 5.034 0.277 
9 3.265 6.371 0.277 
10 4.031 7.865 0.277 
11 4.878 9.517 0.277 
12 5.805 11.326 0.277 
13 6.812 13.292 0.277 
14 7.901 15.416 0.277 
15 9.07 17.697 0.277 




















Ao was obtained from (17) while kọ was obtained from 
(15). The value of deep water depth for each wave period 
is calculated by (18), where kọ is calculated by (15) using 
the weighting coefficient y= 1.211 and 0 = 0.625 as 
depicted in Table (2). 


Deep water depth hoin (18) and in Table 2, only shows the 
water depth limits where the wavenumber or wavelength at 
water depth greater than hAgis a constant value. The 
wavenumber is affected by water depth at a water depth of 
less than ho. 


The value of zo = 0.277, half of the criteria of the Shore 
0 


Protection Manual (1984) which uses a = 0.5. 
0 


V. BREAKER HEIGHT, aa EQUATION 


H ; : f 
z Equation was formulated using the wave amplitude 
b 


function (12). The equation is differentiated on the x- 
horizontal axis for sloping bottom, 
dA B m A BGk dkA 

=O —— Ch — — | — ———— — 


dax 2yo dx 2) 4yo dx 


E was substitude by (5) and“ elaborated, 


aA B L a 2 aa 2 
ax 2ya2dx\" 2) dya\"ax “ax 


f i .. dA 
Breaking occurs in condition =0 
x 
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f _dk kA\ BGk 7 dk 
EE (o-9)- AB 
4yo dx 2 4yo \ dx 


The relation between wavenumber and wave amplitude at 
the breaker point is obtained, 


kpAp = 2 Cq (19) 


Or, 
H, 2 
i, 3n” 


If maximumc,is used,wherec, = 1, therefore: 

H 

—? = 0.212 

Ly 

It was discovered that—is a constant value that applies to 
b 

all wave periods. (19) is substituted into (15) to get the 


breaker wave number equation, 


3 
kp = say o Ry — cf)...(20) 


In this equation, there is no variable of breaker height H, 
or breaker amplitude A, but the formulation used the 


. H ; . ; : 
relation = . Thus, it can be said that in (20), there is an 
b 


effect of breaker wave amplitude, or there is an interaction 
between breaker wavenumber and breaker amplitude. 


The breaker depth equation can be obtained by using the 
wavenumber conservation equation, which is, 


Ap Ao 
ky (hy +2) = ko (ho +) = On 


So, 
Ap 
ky (h + >) = On 
2 
Or, 
= OE An 
ha A ee (21) 


From (19), (20), and (21), it can be said that in the 
breaking equation obtained there is an interaction among 
breaker characteristics. 


VI. DETERMINATION OF THE VALUES OF y 
AND 0 


To get the weighting coefficienty and deep water depth 
coefficient, the breaker height model is calibrated to the 
breaker height from the breaker height index (BHI) and the 
breaker depth model is calibrated to the breaker depth 
from the breaker depth index. 


Breaker height comparison used was the mean value of the 
breaker height of five BHI (Gourlay, 1992; 
Komar&Gaughan,1972; Kraus, 1989; Smith & Kraus, 
1990; Pitikon&Shibayama, 2000). Meanwhile, the breaker 
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depth is compared with the breaker depth from SPM 
(1984). 


Komar and Gaughan (1972) 
1 
Hp = Ho TE 
= 0.56 (=) E 22) 
Larson and Kraus (1989), 
Hp = Ho —0.24 
mb = 0.53 (=) eee (23) 
Smith and Kraus (1990), 
Hp _ Ho —0.30+0.88m 
i = (0.34 + 2.74m) (=) ee (24) 
Gourlay (1992), 
Hp _ Ho -0.28 
ab = 0.478 (=) A ere: (25) 


Pitikon and Shibayama (2000) : 


H, Ho 3 
Fe = (10.02m? — 7.46m? + 1.32m + 0.55) (~) 


0 0 


Hois deep-water wave height, Lo is deep water wavelength 
: ; 2 2 
(calculated using linear wave theory, kọ = ToLo = = m 
0 


is bottom slope and Hy is breaker height. 


Equation for breaker depth index from SPM (1984) is: 
hp 1 Hp 


2 —anorhy, == 71 Fee (27) 
sr) (i) 
z 1.56 
a = 43.75(1 — °°" )b = em 


h, is breaker depth. 

For the deep-water wave height Hy) for the calculation 
input with BHI, the maximum breaker height (17) is used. 
In Tables (3) and (4), the results of calculations and 
comparisons between the results of the model and the 
results of BHI are presented, where the model used 
werey = 1.201and 0 = 0.623. 


Table.3: The Comparison ofBreaker Height H,Model 


























andH), BHI. 
T Ho H, (m) 
(sec) (m) Model BHI 
6 1.451 1.721 1.721 
7 1.975 2.343 2.343 
8 2.58 3.06 3.06 
9 3.265 3.872 3.873 
10 4.031 4.781 4.781 
11 4.878 5.785 5.785 
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12 5.805 6.884 6.885 
13 6.812 8.079 8.08 

14 7.901 9.37 9.371 
15 9.07 10.757 | 10.757 




















In Table (4), the breaker depth index of i: model is 
compatible with z of SPM, and also in accordance with 
Mc Cowan’s equation (1894): 

- = 0.78 ......(28) 


Table.4: The Comparison of breaker depth h Model 
andh, SPM. 





T h, (m) hy 
(sec) | Model SPM model SPM 
6 2.199 2.207 0.783 0.78 
7 2.993 3.003 0.783 0.78 
8 3.909 3.923 0.783 0.78 
9 4.948 4.965 0.783 0.78 
10 6.108 6.129 0.783 0.78 
11 7.391 7.417 0.783 0.78 
12 8.796 8.826 0.783 0.78 
13 10.323 | 10.359 | 0.783 0.78 
14 11.972 | 12.014 | 0.783 0.78 
15 13.743 | 13.791 | 0.783 0.78 
























































There is a relatively small difference in the value of y. In 
deep water, the value of y = 1.211 is obtained by 
calibrating to Michell’s (1893) criteria. While the value of 
y = 1.201is obtained the calibration of breaker height and 
breaker depth. However, this difference is also caused by 
the conditions of the Michell’s criteria (1893) and the 
conditions of the breaker height index. If Michell’s criteria 
are considered more accurate, y = 1.211can be used, 
while if the breaker height index is more accurate, y = 
1.201 can be used. However, with this relatively small 
difference, the mean value of the two can be used, which 
slightly increasing the Michell’s criteria and lowering the 
breaker height. 


VII. THE COMPARISON OF BREAKER LENGTH 


There are a number of breaker length indexes which have 
almost the same basic form. They are: 


a. Miche (1944) 
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Hp _ 2mhp 

T = 0.142 tanh (=) n (29) 
b. Battjes and Jansen (19878) 
H 0.8 2hp 


b 
Hb _ 0.142 tanh = m ) = (30) 


Lp 


Cc. Battjes and Stive (1985) 





H, Ho\ 21h, 
T = 0.14 tanh (0.5 + 0.4 tanh 3372) 


b Lo/ 0.88L, 
eee (31) 

d. Kamphuis (1991) 

Hp _ amy 

72 = 0.127 exp(4m) tanh ( m ) — (32) 

e. Rattanapitikon and Shibayama (2000) 





Hp 2h, 
7 = 0.14 tanh | (-11.21 m? + 5.01 m + 0.91) — 
b b 


...(33) 


In the Breaker Length Index equations, m is the bottom 
slope. In (31), Hp is the deep-water wave height. While Lo 
is the deep-water wave length of the linear wave theory 


2 
dispersion equation of k = > 
The breaker height H, was calculated using the Breaker 
Height Index. The breaker depth is obtained using (28). 


With the input breaker height Hp and breaker depth hy, the 
breaker length can be calculated by (29) — (33). 


Table.5: Breaker Length Index value=” from (29)-(30) 
b 





r m, 
Lp 





(sec) | (29) | (30) | (31) | (32) (33) 

6 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 
0.083 | 0.038 | 0.087 | 0.03 | 0.036 
0.083 | 0.038 | 0.087 | 0.03 | 0.036 
0.083 | 0.038 | 0.087 | 0.03 | 0.036 
10 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 
11 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 
12 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 
13 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 
14 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 
15 | 0.083 | 0.038 | 0.087 | 0.03 | 0.036 














S| oo; N 
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. H . 
The results of the calculation of with a number of 
b 


breaker length index equations vary in value, with a very 
widespread. In addition, it is also much different from the 
critical wave steepness criteria both from Michell (1893) 
and Toffoli et al. (2010) where the criteria of Toffoli et al. 
(2010): 


He = 0.170.....(34) 


Lp 


Therefore, the value of the weighting coefficient y and 
deep water coefficient @ in this study were adjusted to the 
breaker height index and breaker depth index only. 


Comparison of breaker length model with y = 1.201 and 
0 = 0.623, with breaker length (29), is presented in Table 
(6). Where the breaker length of the model is much shorter 


with the breaker length index index{ also much bigger. 
Table.6: The Comparison ofBreaker Length Model 
andMiche’sBreaker Length 

T Lp (m) H, 

Lp 
(sec) | model | Miche | model | Miche 
6 8.44 20.802 | 0.204 0.083 
7 11.488 | 28.314 | 0.204 0.083 
8 15.005 | 36.982 | 0.204 0.083 
9 18.991 | 46.805 | 0.204 0.083 
10 | 23.446 | 57.785 | 0.204 0.083 
11 28.369 | 69.919 | 0.204 0.083 
12 | 33.762 | 83.21 0.204 0.083 
13 39.623 | 97.656 | 0.204 0.083 
14 | 45.954 | 113.258 | 0.204 0.083 
15 | 52.753 | 130.015 | 0.204 0.083 



























































The breaker length index of the model is also much greater 
than the critical wave steepness criteria of Michel (1893) 
and Toffoli et al. (2010). This is due to the wave energy 
compression when the waves move from deeper waters to 
shallower waters making the waves are stronger. 


VIII. CONCLUSION 


The equations for breaking characteristics in this study 
were formulated based on hydrodynamic equations, 
continuity equations, momentum equations, and 
conservation equations attached to the velocity potential, 
which were the energy conservation equation and the 
wavenumber conservation equation. Hence, first, it can be 
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concluded that the obtained equations satisfy the 
conservation laws of fluid dynamics. 


The equations formulated are in line with the previous 
researchers in terms of breaker height and breaker depth. 
Considering that the breaker height is determined by the 
breaker depth and the breaker depth is determined by the 
breaker height and breaker length, it can be concluded that 
the breaker length formulated is appropriate or accurate. 
Thus, in general, it can be concluded that the breaking 
equations formulated can provide good information 
regarding the condition of the breaking wave 
characteristic. 


There are still physical parameters that have not been 
included as a variable in the breaking equations obtained, 
which is the bottom slope. The next research is expected to 
develop a breaking equation involving bottom slope as the 
variable. 
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